PRPF8-mediated dysregulation of hBrr2 helicase disrupts human spliceosome kinetics and 5´-splice-site selection causing tissue-specific defects.

The carboxy-terminus of the spliceosomal protein PRPF8, which regulates the RNA helicase Brr2, is a hotspot for mutations causing retinitis pigmentosa-type 13, with unclear role in human splicing and tissue-specificity mechanism. We used patient induced pluripotent stem cells-derived cells, carrying the heterozygous PRPF8 c.6926 A > C (p.H2309P) mutation to demonstrate retinal-specific endophenotypes comprising photoreceptor loss, apical-basal polarity and ciliary defects. Comprehensive molecular, transcriptomic, and proteomic analyses revealed a role of the PRPF8/Brr2 regulation in 5'-splice site (5'SS) selection by spliceosomes, for which disruption impaired alternative splicing and weak/suboptimal 5'SS selection, and enhanced cryptic splicing, predominantly in ciliary and retinal-specific transcripts. Altered splicing efficiency, nuclear speckles organisation, and PRPF8 interaction with U6 snRNA, caused accumulation of active spliceosomes and poly(A)+ mRNAs in unique splicing clusters located at the nuclear periphery of photoreceptors. Collectively these elucidate the role of PRPF8/Brr2 regulatory mechanisms in splicing and the molecular basis of retinal disease, informing therapeutic approaches.

A novel homozygous splice site variant in ARL2BP causes a syndromic autosomal recessive rod-cone dystrophy with situs inversus, asthenozoospermia, unilateral renal agenesis and microcysts.

BACKGROUND: This report presents a clinical case of syndromic rod-cone dystrophy due to a splice site variant in the ARL2BP gene causing situs inversus, asthenozoospermia, unilateral renal agenesis and microcysts. The presence of renal agenesis and cryptorchidism expands the clinical manifestations due to ARL2BP variants. The detailed, long-term follow-up contributes valuable insights into disease progression, aiding clinical diagnosis and patient management. CASE PRESENTATION: The male patient complained of photophobia as the first symptom when he was 20 years old followed by nyctalopia, loss of central visual acuity and peripheral visual field ten years later. Genetic analysis identified a likely pathogenic homozygous variant (c.294-1G > C) involving the splicing acceptor site of intron 4. Reported symptoms together with full-field stimulus threshold testing, electroretinogram and advanced multimodal imaging allowed us to recognize the typical characteristics of a mixed retinal dystrophy. Despite the end-stage retinal disease, this patient still retained a useful residual vision at 63 years and had a slow disease progression during the last 5 years of evaluation. DISCUSSION AND CONCLUSIONS: Our findings underscore the variable clinical presentation of ARL2BP variants, emphasizing the importance of a nuanced approach in diagnosing and managing patients. The presence of renal cysts warrants consideration of a differential diagnosis, particularly with Senior-Loken (SLS), Bardet-Biedl (BBS) and Joubert syndromes (JS) but also with Short Rib Thoracic Dysplasia 9, highlighting the need for careful phenotypic evaluation in these cases.

Structural basis of U12-type intron engagement by the fully assembled human minor spliceosome.

The minor spliceosome, which is responsible for the splicing of U12-type introns, comprises five small nuclear RNAs (snRNAs), of which only one is shared with the major spliceosome. In this work, we report the 3.3-angstrom cryo-electron microscopy structure of the fully assembled human minor spliceosome pre-B complex. The atomic model includes U11 small nuclear ribonucleoprotein (snRNP), U12 snRNP, and U4atac/U6atac.U5 tri-snRNP. U11 snRNA is recognized by five U11-specific proteins (20K, 25K, 35K, 48K, and 59K) and the heptameric Sm ring. The 3' half of the 5'-splice site forms a duplex with U11 snRNA; the 5' half is recognized by U11-35K, U11-48K, and U11 snRNA. Two proteins, CENATAC and DIM2/TXNL4B, specifically associate with the minor tri-snRNP. A structural analysis uncovered how two conformationally similar tri-snRNPs are differentiated by the minor and major prespliceosomes for assembly.

A splice acceptor variant in RGS6 associated with intellectual disability, microcephaly, and cataracts disproportionately promotes expression of a subset of RGS6 isoforms.

Intellectual disability (ID) is associated with an increased risk of developing psychiatric disorders, suggesting a common underlying genetic factor. Importantly, altered signaling and/or expression of regulator of G protein signaling 6 (RGS6) is associated with ID and numerous psychiatric disorders. RGS6 is highly conserved and undergoes complex alternative mRNA splicing producing ~36 protein isoforms with high sequence similarity historically necessitating a global approach in functional studies. However, our recent analysis in mice revealed RGS6 is most highly expressed in CNS with RGS6L(+GGL) isoforms predominating. A previously reported genetic variant in intron 17 of RGS6 (c.1369-1G>C), associated with ID, may provide further clues into RGS6L(+GGL) isoform functional delineation. This variant was predicted to alter a highly conserved canonical 3' acceptor site creating an alternative branch point within exon 18 (included in a subset of RGS6L(+GGL) transcripts) and a frameshift forming an early stop codon. We previously identified this alternative splice site and demonstrated its use generates RGS6Lζ(+GGL) isoforms. Here, we show that the c.1369-1G>C variant disrupts the canonical, preferred (>90%) intron 17 splice site and leads to the exclusive use of the alternate exon 18 splice site, inducing disproportionate expression of a subset of isoforms, particularly RGS6Lζ(+GGL). Furthermore, RGS6 global knockout mice do not exhibit ID. Thus, ID caused by the c.1369-1G>C variant likely results from altered RGS6 isoform expression, rather than RGS6 isoform loss. In summary, these studies highlight the importance of proper RGS6 splicing and identify a previously unrecognized role of G protein signaling in ID.

A noncanonical splicing variant c.875-5 T > G in von Willebrand factor causes in-frame exon skipping and type 2A von Willebrand disease.

INTRODUCTION: A novel variant involving noncanonical splicing acceptor site (c.875-5 T > G) in propeptide coding region of von Willebrand factor (VWF) was identified in a patient with type 2A von Willebrand disease (VWD), who co-inherited with a null variant (p.Tyr271*) and presented characteristic discrepancy of plasma level of VWF antigen and activity, and a selective reduction of both intermediate-molecular-weight (IMWMs) and high-molecular-weight VWF multimers (HMWMs). MATERIALS AND METHODS: VWF mRNA transcripts obtained from peripheral leukocytes and platelets of the patients were investigated to analyze the consequence of c.875-5 T > G on splicing. The impact of the variant on expression and multimer assembly was further analyzed by in vitro expression studies in AtT-20 cells. The intracellular processing of VWF mutant and the Weibel-Palade bodies (WPBs) formation was evaluated by immunofluorescence staining and electron microscopy. RESULTS: The mRNA transcript analysis revealed that c.875-5 T > G variant led to exon 8 skipping and an in-frame deletion of 41 amino acids in the D1 domain of VWF (p.Ser292_Glu333delinsLys), yielding a truncated propeptide. Consistent with the patient's laboratory manifestations, the AtT-20 cells transfected with mutant secreted less VWF, with the VWF antigen level in conditioned medium 47 % of wild-type. A slight retention in the endoplasmic reticulum was observed for the mutant. Almost complete loss of IMWMs and HMWMs in the medium and impaired WPBs formation in the cell, indicating truncated VWF propeptide lost its chaperon-like function for VWF multimerization and tubular storage. CONCLUSIONS: The VWF splicing site variant (c.875-5 T > G) causes propeptide truncation, severely compromising VWF multimer assembly and tubular storage.

Cryo-EM analyses of dimerized spliceosomes provide new insights into the functions of B complex proteins.

The B complex is a key intermediate stage of spliceosome assembly. To improve the structural resolution of monomeric, human spliceosomal B (hB) complexes and thereby generate a more comprehensive hB molecular model, we determined the cryo-EM structure of B complex dimers formed in the presence of ATP γ S. The enhanced resolution of these complexes allows a finer molecular dissection of how the 5' splice site (5'ss) is recognized in hB, and new insights into molecular interactions of FBP21, SNU23 and PRP38 with the U6/5'ss helix and with each other. It also reveals that SMU1 and RED are present as a heterotetrameric complex and are located at the interface of the B dimer protomers. We further show that MFAP1 and UBL5 form a 5' exon binding channel in hB, and elucidate the molecular contacts stabilizing the 5' exon at this stage. Our studies thus yield more accurate models of protein and RNA components of hB complexes. They further allow the localization of additional proteins and protein domains (such as SF3B6, BUD31 and TCERG1) whose position was not previously known, thereby uncovering new functions for B-specific and other hB proteins during pre-mRNA splicing.

Specificity, synergy, and mechanisms of splice-modifying drugs.

Drugs that target pre-mRNA splicing hold great therapeutic potential, but the quantitative understanding of how these drugs work is limited. Here we introduce mechanistically interpretable quantitative models for the sequence-specific and concentration-dependent behavior of splice-modifying drugs. Using massively parallel splicing assays, RNA-seq experiments, and precision dose-response curves, we obtain quantitative models for two small-molecule drugs, risdiplam and branaplam, developed for treating spinal muscular atrophy. The results quantitatively characterize the specificities of risdiplam and branaplam for 5' splice site sequences, suggest that branaplam recognizes 5' splice sites via two distinct interaction modes, and contradict the prevailing two-site hypothesis for risdiplam activity at SMN2 exon 7. The results also show that anomalous single-drug cooperativity, as well as multi-drug synergy, are widespread among small-molecule drugs and antisense-oligonucleotide drugs that promote exon inclusion. Our quantitative models thus clarify the mechanisms of existing treatments and provide a basis for the rational development of new therapies.

Large-scale analysis reveals splicing biomarkers for tuberculosis progression and prognosis.

BACKGROUND: Emerging evidence suggests that aberrant alternative splicing (AS) may play an important role in tuberculosis (TB). However, current knowledge regarding the value of AS in TB progression and prognosis remains unclear. METHOD: Public RNA-seq datasets related to TB progression and prognosis were searched and AS analyses were conducted based on SUPPA2. Percent spliced in (PSI) was used for quantifying AS events and multiple machine learning (ML) methods were employed to construct predictive models. Area under curve (AUC), sensitivity and specificity were calculated to evaluate the model performance. RESULTS: A total of 1587 samples from 7 datasets were included. Among 923 TB-progression related differential AS events (DASEs), 3 events (GET1-skipping exon (SE), TPD52-alternative first exons (AF) and TIMM10-alternative 5' splice site (A5)) were selected as candidate biomarkers; however, their predictive performance was limited. For TB prognosis, 5 events (PHF23-AF, KIF1B-SE, MACROD2-alternative 3' splice site (A3), CD55-retained intron (RI) and GALNT11-AF) were selected as candidates from the 1282 DASEs. Six ML methods were used to integrate these 5 events and XGBoost outperformed than others. AUC, sensitivity and specificity of XGBoost model were 0.875, 81.1% and 83.5% in training set, while they were 0.805, 68.4% and 73.2% in test set. CONCLUSION: GET1-SE, TPD52-AF and TIMM10-A5 showed limited role in predicting TB progression, while PHF23-AF, KIF1B-SE, MACROD2-A3, CD55-RI and GALNT11-AF could well predict TB prognosis and work as candidate biomarkers. This work preliminarily explored the value of AS in predicting TB progression and prognosis and offered potential targets for further research.

Spliceosomal helicases DDX41/SACY-1 and PRP22/MOG-5 both contribute to proofreading against proximal 3' splice site usage.

RNA helicases drive necessary rearrangements and ensure fidelity during the pre-mRNA splicing cycle. DEAD-box helicase DDX41 has been linked to human disease and has recently been shown to interact with DEAH-box helicase PRP22 in the spliceosomal C* complex, yet its function in splicing remains unknown. Depletion of DDX41 homolog SACY-1 from somatic cells has been previously shown to lead to changes in alternative 3' splice site (3'ss) usage. Here, we show by transcriptomic analysis of published and novel data sets that SACY-1 perturbation causes a previously unreported pattern in alternative 3' splicing in introns with pairs of 3' splice sites ≤18 nt away from each other. We find that both SACY-1 depletion and the allele sacy-1(G533R) lead to a striking unidirectional increase in the usage of the proximal (upstream) 3'ss. We previously discovered a similar alternative splicing pattern between germline tissue and somatic tissue, in which there is a unidirectional increase in proximal 3'ss usage in the germline for ∼200 events; many of the somatic SACY-1 alternative 3' splicing events overlap with these developmentally regulated events. We generated targeted mutant alleles of the Caenorhabditis elegans homolog of PRP22, mog-5, in the region of MOG-5 that is predicted to interact with SACY-1 based on the human C* structure. These viable alleles, and a mimic of the myelodysplastic syndrome-associated allele DDX41(R525H), all promote the usage of proximal alternative adjacent 3' splice sites. We show that PRP22/MOG-5 and DDX41/SACY-1 have overlapping roles in proofreading the 3'ss.

A hypomorphic variant of choroideremia is associated with a novel intronic mutation that leads to exon skipping.

INTRODUCTION: Molecular confirmation of pathogenic sequence variants in the CHM gene is required prior to enrolment in retinal gene therapy clinical trials for choroideremia. Individuals with mild choroideremia have been reported. The molecular basis of genotype-phenotype associations is of clinical relevance since it may impact on selection for retinal gene therapy. METHODS AND MATERIALS: Genetic testing and RNA analysis were undertaken in a patient with mild choroideremia to confirm the pathogenicity of a novel intronic variant in CHM and to explore the mechanism underlying the mild clinical phenotype. RESULTS: A 42-year-old male presented with visual field loss. Fundoscopy and autofluorescence imaging demonstrated mild choroideremia for his age. Genetic analysis revealed a variant at a splice acceptor site in the CHM gene (c.1350-3C > G). RNA analysis demonstrated two out-of-frame transcripts, suggesting pathogenicity, without any detectable wildtype transcripts. One of the two out-of-frame transcripts is present in very low levels in healthy controls. DISCUSSION: Mild choroideremia may result from +3 or -3 splice site variants in CHM. It is presumed that the resulting mRNA transcripts may be partly functional, thereby preventing the development of the null phenotype. Choroideremia patients with such variants may present challenges for gene therapy since there may be residual transcript activity which could result in long-lasting visual function which is atypical for this disease.

Pleiotropic effects of different exonic nucleotide changes at the same position contribute to hemophilia B phenotypic variation.

BACKGROUND: The disease-causing effects of genetic variations often depend on their location within a gene. Exonic changes generally lead to alterations in protein production, secretion, activity, or clearance. However, owing to the overlap between proteins and splicing codes, missense variants can also affect messenger RNA splicing, thus adding a layer of complexity and influencing disease phenotypes. OBJECTIVES: To extensively characterize a panel of 13 exonic variants in the F9 gene occurring at 6 different factor IX positions and associated with varying severities of hemophilia B (HB). METHODS: Computational predictions, splicing analysis, and recombinant factor IX assays were exploited to characterize F9 variants. RESULTS: We demonstrated that 5 (38%) of 13 selected F9 exonic variants have pleiotropic effects. Although bioinformatic approaches accurately classified effects, extensive experimental assays were required to elucidate and deepen the molecular mechanisms underlying the pleiotropic effects. Importantly, their characterization was instrumental in developing tailored RNA therapeutics based on engineered U7 small nuclear RNA to mask cryptic splice sites and compensatory U1 small nuclear RNA to enhance exon definition. CONCLUSION: Overall, albeit a multitool bioinformatic approach suggested the molecular effects of multiple HB variants, the deep investigation of molecular mechanisms revealed insights into the HB phenotype-genotype relationship, enabling accurate classification of HB variants. Importantly, knowledge of molecular mechanisms allowed the development of tailored RNA therapeutics, which can also be translated to other genetic diseases.

CAG repeat expansions create splicing acceptor sites and produce aberrant repeat-containing RNAs.

Expansions of CAG trinucleotide repeats cause several rare neurodegenerative diseases. The disease-causing repeats are translated in multiple reading frames and without an identifiable initiation codon. The molecular mechanism of this repeat-associated non-AUG (RAN) translation is not known. We find that expanded CAG repeats create new splice acceptor sites. Splicing of proximal donors to the repeats produces unexpected repeat-containing transcripts. Upon splicing, depending on the sequences surrounding the donor, CAG repeats may become embedded in AUG-initiated open reading frames. Canonical AUG-initiated translation of these aberrant RNAs may account for proteins that have been attributed to RAN translation. Disruption of the relevant splice donors or the in-frame AUG initiation codons is sufficient to abrogate RAN translation. Our findings provide a molecular explanation for the abnormal translation products observed in CAG trinucleotide repeat expansion disorders and add to the repertoire of mechanisms by which repeat expansion mutations disrupt cellular functions.

The unusual gene architecture of polyubiquitin is created by dual-specific splice sites.

BACKGROUND: The removal of introns occurs through the splicing of a 5' splice site (5'ss) with a 3' splice site (3'ss). These two elements are recognized by distinct components of the spliceosome. However, introns in higher eukaryotes contain many matches to the 5' and 3' splice-site motifs that are presumed not to be used. RESULTS: Here, we find that many of these sites can be used. We also find occurrences of the AGGT motif that can function as either a 5'ss or a 3'ss-previously referred to as dual-specific splice sites (DSSs)-within introns. Analysis of the Sequence Read Archive reveals a 3.1-fold enrichment of DSSs relative to expectation, implying synergy between the ability to function as a 5'ss and 3'ss. Despite this suggested mechanistic advantage, DSSs are 2.7- and 4.7-fold underrepresented in annotated 5' and 3' splice sites. A curious exception is the polyubiquitin gene UBC, which contains a tandem array of DSSs that precisely delimit the boundary of each ubiquitin monomer. The resulting isoforms splice stochastically to include a variable number of ubiquitin monomers. We found no evidence of tissue-specific or feedback regulation but note the 8.4-fold enrichment of DSS-spliced introns in tandem repeat genes suggests a driving role in the evolution of genes like UBC. CONCLUSIONS: We find an excess of unannotated splice sites and the utilization of DSSs in tandem repeats supports the role of splicing in gene evolution. These findings enhance our understanding of the diverse and complex nature of the splicing process.

Structural basis of branching during RNA splicing.

Branching is a critical step in RNA splicing that is essential for 5' splice site selection. Recent spliceosome structures have led to competing models for the recognition of the invariant adenosine at the branch point. However, there are no structures of any splicing complex with the adenosine nucleophile docked in the active site and positioned to attack the 5' splice site. Thus we lack a mechanistic understanding of adenosine selection and splice site recognition during RNA splicing. Here we present a cryo-electron microscopy structure of a group II intron that reveals that active site dynamics are coupled to the formation of a base triple within the branch-site helix that positions the 2'-OH of the adenosine for nucleophilic attack on the 5' scissile phosphate. This structure, complemented with biochemistry and comparative analyses to splicing complexes, supports a base triple model of adenosine recognition for branching within group II introns and the evolutionarily related spliceosome.

MAJIQlopedia: an encyclopedia of RNA splicing variations in human tissues and cancer.

Quantification of RNA splicing variations based on RNA-Sequencing can reveal tissue- and disease-specific splicing patterns. To study such splicing variations, we introduce MAJIQlopedia, an encyclopedia of splicing variations that encompasses 86 human tissues and 41 cancer datasets. MAJIQlopedia reports annotated and unannotated splicing events for a total of 486 175 alternative splice junctions in normal tissues and 338 317 alternative splice junctions in cancer. This database, available at https://majiq.biociphers.org/majiqlopedia/, includes a user-friendly interface that provides graphical representations of junction usage quantification for each junction across all tissue or cancer types. To demonstrate case usage of MAJIQlopedia, we review splicing variations in genes WT1, MAPT and BIN1, which all have known tissue or cancer-specific splicing variations. We also use MAJIQlopedia to highlight novel splicing variations in FDX1 and MEGF9 in normal tissues, and we uncover a novel exon inclusion event in RPS6KA6 that only occurs in two cancer types. Users can download the database, request the addition of data to the webtool, or install a MAJIQlopedia server to integrate proprietary data. MAJIQlopedia can serve as a reference database for researchers seeking to understand what splicing variations exist in genes of interest, and those looking to understand tissue- or cancer-specific splice isoform usage.

Functional analysis of the CTNS gene exonic variants predicted to affect splicing.

Cystinosis is a severe, monogenic systemic disease caused by variants in CTNS gene. Currently, there is growing evidence that exonic variants in many diseases can affect pre-mRNA splicing. The impact of CTNS gene exonic variants on splicing regulation may be underestimated due to the lack of routine studies at the RNA level. Here, we analyzed 59 exonic variants in the CTNS gene using bioinformatics tools and identified candidate variants that may induce splicing alterations by minigene assays. We identified six exonic variants that induce splicing alterations by disrupting the ratio of exonic splicing enhancers/exonic splicing silencers (ESEs/ESSs) or by interfering with the recognition of classical splice sites, or both. Our results help in the correct molecular characterization of variants in cystinosis and inform emerging therapies. Furthermore, our work suggests that the combination of in silico and in vitro assays facilitates to assess the effects of DNA variants driving rare genetic diseases on splicing regulation and will enhance the clinical utility of variant functional annotation.

Estimating the proportion of nonsense variants undergoing the newly described phenomenon of manufactured splice rescue.

A recent report described a nonsense variant simultaneously creating a donor splice site, resulting in a truncated but functional protein. To explore the generalizability of this unique mechanism, we annotated >115,000 nonsense variants using SpliceAI. Between 0.61% (donor gain delta score >0.8, for high precision) and 2.57% (>0.2, for high sensitivity) of nonsense variants were predicted to create new donor splice sites at or upstream of the stop codon. These variants were less likely than other nonsense variants in the same genes to be classified as pathogenic/likely pathogenic in ClinVar (p < 0.001). Up to 1 in 175 nonsense variants were predicted to result in small in-frame deletions and loss-of-function evasion through this "manufactured splice rescue" mechanism. We urge caution when interpreting nonsense variants where manufactured splice rescue is a strong possibility and correlation with phenotype is challenging, as will often be the case with secondary findings and newborn genomic screening programs.

Exonic splicing code and coordination of divalent metals in proteins.

Exonic sequences contain both protein-coding and RNA splicing information but the interplay of the protein and splicing code is complex and poorly understood. Here, we have studied traditional and auxiliary splicing codes of human exons that encode residues coordinating two essential divalent metals at the opposite ends of the Irving-Williams series, a universal order of relative stabilities of metal-organic complexes. We show that exons encoding Zn2+-coordinating amino acids are supported much less by the auxiliary splicing motifs than exons coordinating Ca2+. The handicap of the former is compensated by stronger splice sites and uridine-richer polypyrimidine tracts, except for position -3 relative to 3' splice junctions. However, both Ca2+ and Zn2+ exons exhibit close-to-constitutive splicing in multiple tissues, consistent with their critical importance for metalloprotein function and a relatively small fraction of expendable, alternatively spliced exons. These results indicate that constraints imposed by metal coordination spheres on RNA splicing have been efficiently overcome by the plasticity of exon-intron architecture to ensure adequate metalloprotein expression.

Systematic Minigene-Based Splicing Analysis and Tentative Clinical Classification of 52 CHEK2 Splice-Site Variants.

BACKGROUND: Disrupted pre-mRNA splicing is a frequent deleterious mechanism in hereditary cancer. We aimed to functionally analyze candidate spliceogenic variants of the breast cancer susceptibility gene CHEK2 by splicing reporter minigenes. METHODS: A total of 128 CHEK2 splice-site variants identified in the Breast Cancer After Diagnostic Gene Sequencing (BRIDGES) project (https://cordis.europa.eu/project/id/634935) were analyzed with MaxEntScan and subsetted to 52 variants predicted to impact splicing. Three CHEK2 minigenes, which span all 15 exons, were constructed and validated. The 52 selected variants were then genetically engineered into the minigenes and assayed in MCF-7 (human breast adenocarcinoma) cells. RESULTS: Of 52 variants, 46 (88.5%) impaired splicing. Some of them led to complex splicing patterns with up to 11 different transcripts. Thirty-four variants induced splicing anomalies without any trace or negligible amounts of the full-length transcript. A total of 89 different transcripts were annotated, which derived from different events: single- or multi-exon skipping, alternative site-usage, mutually exclusive exon inclusion, intron retention or combinations of the abovementioned events. Fifty-nine transcripts were predicted to introduce premature termination codons, 7 kept the original open-reading frame, 5 removed the translation start codon, 6 affected the 5'UTR (Untranslated Region), and 2 included missense variations. Analysis of variant c.684-2A > G revealed the activation of a non-canonical TG-acceptor site and exon 6 sequences critical for its recognition. CONCLUSIONS: Incorporation of minigene read-outs into an ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based classification scheme allowed us to classify 32 CHEK2 variants (27 pathogenic/likely pathogenic and 5 likely benign). However, 20 variants (38%) remained of uncertain significance, reflecting in part the complex splicing patterns of this gene.

Global analysis of binding sites of U2AF1 and ZRSR2 reveals RNA elements required for mutually exclusive splicing by the U2- and U12-type spliceosome.

Recurring mutations in genes encoding 3' splice-site recognition proteins, U2AF1 and ZRSR2 are associated with human cancers. Here, we determined binding sites of the proteins to reveal that U2-type and U12-type splice sites are recognized by U2AF1 and ZRSR2, respectively. However, some sites are spliced by both the U2-type and U12-type spliceosomes, indicating that well-conserved consensus motifs in some U12-type introns could be recognized by the U2-type spliceosome. Nucleotides flanking splice sites of U12-type introns are different from those flanking U2-type introns. Remarkably, the AG dinucleotide at the positions -1 and -2 of 5' splice sites of U12-type introns with GT-AG termini is not present. AG next to 5' splice site introduced by a single nucleotide substitution at the -2 position could convert a U12-type splice site to a U2-type site. The class switch of introns by a single mutation and the bias against G at the -1 position of U12-type 5' splice site support the notion that the identities of nucleotides in exonic regions adjacent to splice sites are fine-tuned to avoid recognition by the U2-type spliceosome. These findings may shed light on the mechanism of selectivity in U12-type intron splicing and the mutations that affect splicing.